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(57) Abstract 



A device and method for precoding data signals for pulse code modulation (PCM) transmission including a transmitter (52) for 
transmitting a sequence of analog levels over analog channel to a quantization device, wherein the analog channel modifies the transmitted 
analog levels, the transmitter (52) comprising: a mapping device (150) for mapping data bits to be transmitted to a sequence of equivalence 
classes, wherein each equivalence class contains one or more constellation points; and a constellation point selector (152) interconnected 
to the mapping device (150) which selects a constellation point in each equivalence class to represent the data bits to be transmitted and 
which transmits an analog level that produces the selected constellation point at an input to the quantization device. 
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DEVICE AND METHOD FOR PRECODING DATA SIGNALS 
FOR PCM TRANSMISSION 

RELATED APPLICATIONS 

This application is a continuation-in-part of U.S. Apl. Ser. No. 08/747,840, filed 
November 13, 1996, which is hereby incorporated by reference in its entirety. 

FIELD OF INVENTION 

This invention relates to a device and method for precoding data signals for 
pulse code modulation (PCM) transmission. 

BACKGROUND OF INVENTION 

Conventional modems, such as V.34 modems, treat the public switched 
telephone network (PSTN) as a pure analog channel even though the signals are 
digitized throughout most of the network. In contrast, pulse code modulation (PCM) 
modems take advantage of the fact that most of the network is digital and that 
typically central site modems, such as those of internet service providers and on-line 
services, are connected to the PSTN via digital connections (e.g., T1 in the United 
States and E1 in Europe). First generation PCM modems transmit data in PCM 
mode downstream only (i.e., from a central site digital modem to an analog end user 
modem) and transmit in analog mode, e.g. V.34 mode, upstream (i.e., from the end 
user modem to the central site modem). Future generation PCM modems will also 
transmit data upstream in PCM mode. 

With PCM downstream, the central site PCM modem transmits over a digital 
network eight bit digital words (octets) corresponding to different central office codec 
output levels. At the end user's central office, the octets are converted to analog 
levels which are transmitted over an analog loop. The end user's PCM modem then 
converts the analog levels, viewed as a pulse code amplitude modulated (PAM) 
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signal, into equalized digital levels. The equalized digital levels are ideally mapped 
back into the originally transmitted octets and the data the octets represent. 

With PCM upstream, the end user PCM modem transmits analog levels over 
the analog loop corresponding to data to be transmitted. The analog levels are 
5 modified by the channel characteristics of the analog loop and the modified levels 
are quantized to form octets by a codec in the end user's central office. The codec 
transmits the octets to the PCM central site modem over the digital network. The 
PCM central site modem determines from the octets the transmitted levels and from 
the levels the data transmitted by the end user PCM modem is recovered. 

10 A difficulty that exists with upstream PCM transmission is that the levels 

transmitted by the end user PCM modem are modified by the analog loop. Since 
these modified levels are the levels that are quantized to form octets by the codec, 
and not the levels that are actually transmitted, it can be difficult for the central site 
modem to accurately determine from the octets the data being transmitted by the 

15 end user PCM modem. This difficulty is compounded by the fact that there is a 
channel null in the analog loop, quantization noise introduced by the codec in the 
end user's central office and downstream PCM echo, which make it more difficult for 
the central site PCM modem to accurately recover the data transmitted. 

Therefore, a need exists for a device and method for precoding data signals 

20 for PCM transmission such that the analog levels that are transmitted by the end 
user PCM modem accurately produce predetermined analog levels (constellation 
points) at the input to the codec in the end user's central office, which analog levels 
(constellation points) correspond to the data to be transmitted by the end user PCM 
modem. Moreover, there is a need for a device, system and method for precoding 

25 data signals for PCM transmission which limits the transmit power and combats a 
channel null introduced by the analog loop and quantization noise introduced by the 
codec in the end user's central office. 

BRIEF DESCRIPTION OF THE DRAWINGS 

30 FIG. 1 is a simplified block diagram of a typical telephone company central 

office; 
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FIG. 2 is plot of the frequency spectrum of the yk signals output from the |n-law 
to linear converter of FIG, 1 and the spectral shape of the low pass filter of FIG. 1 ; 

FIG. 3 is a plot of a portion of two frequency spectrums each having a null at 
DC, wherein one spectrum falls off to zero very abruptly at DC and the other 
5 spectrum falls off more gradually; 

FIG. 4 is a diagrammatic representation of a portion of a typical ju,-law 
constellation; 

FIG. 5 is a block diagram of a modem data connection over the telephone 
system including a transmitter for spectrally shaping signals according to this 
10 invention; 

FIG. 6 is a block diagram of the encoder of FIG. 6 used specifically for 
creating a DC null in said analog signals over an analog loop of the PSTN; 

FIG. 7 is a block diagram of the encoder of FIG. 6 which may be used 
generally for modifying, as desired, the frequency spectrum of the signals output 
15 from the analog loop to the end user; 

FIG. 8 is a block diagram of a typical analog PCM modem to digital PCM 
modem communication system; 

FIG. 9 is a more detailed block diagram depicting PCM upstream transmission 
according to this invention; 
20 FIG. 1 0 is an equivalent discrete time block diagram of the block diagram of 

FIG. 9; 

FIG. 1 1 is the equivalent discrete time block diagram of the block diagram of 
FIG. 9 with the analog modem sampling rate twice that of the CO sampling rate; 

FIG. 12 is an example of a transmit constellation having equivalence classes 
25 according to this invention; 

FIG. 13 is a more detailed block diagram of the analog PCM modem 
transmitter of FIG. 10 according to this invention; 

FIG. 14A is another example of a transmit constellation having equivalence 
classes according to this invention; 
30 FIG. 14B is yet another example of a transmit constellation having 

equivalence classes according to this invention; 
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FIG. 15 is a block diagram of a typical analog PCM modem to analog PCM 
modem communication system; 

FIG. 16 is a more detailed block diagram depicting PCM transmission with the 
PCM modem communication system of FIG. 15; and 
5 FIG. 17 is an equivalent discrete time block diagram of the block diagram of 

FIG. 16. 

DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT 

There is first described below a technique for PCM downstream spectral 
10 shaping or precoding of data signals. Then, there is described a precoding 

technique for PCM upstream transmission of data signals. Finally, it is described 
how the PCM upstream precoding technique according to this invention may be 
generalized for use in a PCM communication system interconnecting two analog 
PCM modems, as opposed to the typical analog PCM modem and digital PCM 
15 modem interconnection. 

PCM Downstream Spectral Shapinq/Precodinq 

FIGS. 1 and 2 illustrate the presence of energy near DC in the signals 
transmitted to a remote user's modem over an analog loop. There is shown in FIG. 1 

20 a portion of a typical telephone central office 10 on a PSTN which receives at input 
12 |n-law octets transmitted from a modem (transmitting modem, not shown) directly 
attached to the digital portion of the telephone system, such as the one described in 
the co-pending applications referred to above which directly encodes the digital data 
into octets for transmission. These octets are converted by a D/A converter, also 

25 known as a jj,-law to linear converter 14, to a sequence of voltage levels, yk, each 
level being one of 255 p,-law levels. The levels are output over line 16 to a LPF 18 
which outputs over analog loop 20 towards the remote modem's receiver a filtered 
analog signal s(t) which is an analog representation of the levels. The analog signal 
is demodulated and decoded by the receiving modem which outputs a digital 

30 bitstream which is an estimate of the originally transmitted data. 
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The sequence of levels yk on line 16 from |i-law to linear converter 14 has a 
flat frequency response 22, FIG. 2. The spectral shape 24 of LPF 18 contains a 
significant amount of energy near DC (f=0) as illustrated at point 26. Since the 
sequence yk has a flat frequency response, the spectrum of the signal s(t) output by 
5 filter 18 has the same spectral shape 24 as the filter 18 and therefore the signal s(t) 
also contains a significant amount of energy near DC. As described above, this 
energy near DC tends to saturate the transformers on the system which produces 
unwanted non-linear distortion in the signal s(t) transmitted towards the receiving 
modem. 

10 In some applications this distortion must be reduced. This can be 

accomplished by reducing the signal energy near DC in the transmitted signal. Such 
a DC null 28 is depicted in FIG. 3. As is known in the state-of-the-art, in order to 
create this spectral null at DC in the transmitted signal, the running digital sum (RDS) 
of the transmitted levels yk (namely, the algebraic sum of all previously transmitted 

15 levels) must be kept close to zero. The shape of the spectrum around the DC null 28 
can vary from a relatively shallow sloped spectrum 30 to a spectrum 32 which falls 
off very abruptly at DC, The sharpness of the null depends on how tightly the RDS is 
controlled. 

The present invention accordingly encodes the digital data being transmitted 
20 into p,-law octets in a manner that maintains the RDS near zero to create the desired 
spectral null at DC thereby reducing the non-linear distortion caused by transformer 
saturation. 

To illustrate the method of creating a spectral null, we consider an example of 
transmitting 6 bits with every symbol yk. It will be apparent to those skilled in the art 

25 that the invention can be used for transmitting any other number of bits per symbol, 
or when the number of bits per symbol transmitted varies from symbol to symbol. In 
a system without a spectral null, one first selects a subset of 64 levels from the 
available 255 p.-law levels such that a minimum distance dmin between levels is 
maintained. These 64 levels are symmetric in the sense that for every positive level 

30 there is a negative level of the same magnitude. For example, one can achieve a 
dmin of 32 for an average energy well under -12 dBmO, the regulatory limit. 
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A partial representation of all 255 p,-law levels 34 (128 positive and 127 
negative) is shown in FIG. 4. These levels follow a logarithmic law, with the 64 levels 
closest to the origin being uniformly spaced between -63 and 63 with a spacing of 2. 
The next positive and negative segments start at +/- 66 and they each contain 1 6 
5 points spaced by 4. The scale continues with segments of 16 points, each with a 

spacing of the form 2 n separated from the previous segment by a spacing of .75 * 

2 n . The final segments extend between +/-21 12 and +/-4032 with a spacing of 128. 
The set 35 is the set of 64 levels selected from these 255 levels to represent each 

6 

combination of six bits, i.e. 2 = 64. 

10 In the transmitter, incoming bits are collected in groups of 6, and then mapped 

into p,-law octets, which represent the desired level. In the central office, the p,-law 
octets are converted into levels, and the resulting levels are then transmitted. In the 
receiver, an equalizer compensates for the distortion introduced by the LPF and the 
local loop, and then a decision device estimates the transmitted level, by selecting 

15 the level that is closest to the received point. 

In order to achieve spectral shaping in the above example, additional levels 
are also used, but the minimum distance between levels is still kept at 32. For 
example, consider the case where 92 levels are used. First, these 92 levels are 
divided into equivalence classes. There are a number of different ways for 

20 generating these equivalence classes. One particularly useful way is described 
here: we label the levels by integers 0 through 91 , for example by assigning the 
label 0 to the smallest (most negative) level, the label 1 to the next smallest level, 
and so on. Then, we define 64 "equivalence classes" by grouping together levels 
whose labels differ exactly by 64. Such grouping leads to 36 equivalence classes 

25 with only one level corresponding to one of 36 innermost levels of smallest 

magnitude, and 28 equivalence classes with two levels whose labels differ by 64. 
Other methods for generating the equivalence classes may be used. Each possible 
combination of 6 bits to be transmitted is then represented by an equivalence class. 
For example, the bit combination 000000 may correspond to the first 

30 equivalence class which consists of two levels 
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each being represented by a different octet. Note that it is not necessary to use the 
full dynamic range of the D/A converter. The technique can work with any number of 
levels, as long as more than 64 levels are used. Of course, the more levels used, 
the better the desired spectral shape can be achieved. Our experiments indicate 
5 that very few additional levels need to be considered for generating a DC null with a 
relatively sharp notch. 

In the above example, since each combination of six information bits is 
represented by an equivalence class and often there is more than one level in an 
equivalence class, the information bits must be mapped into one of the levels in a 

10 selected equivalence class before an octet representing that level is transmitted. 
This function is described below with regard to FIGS. 5-7. 

Transmitter 52, FIG. 5, receives from a digital data source, such as a 
computer, a bitstream of digital data and with bit collector 54 divides the bits into 
groups of six, for example. Each six-bit group is provided to encoder 56 which 

15 selects the equivalence classes from which the desired levels to achieve the spectral 
null at DC will be selected. The octets which represent the selected levels are output 
from encoder 56, transmitted over digital circuit-switched telephone network 58 and 
arrive at the remote user's central office 60. At central office 60, the octets are 
converted by p,-law to linear converter 62 to the levels, yk , which pass through LPF 

20 64 and are output over local analog loop 65 as a signal s(t) having a spectral null at 
DC. In receiver 66, the signal s(t) is sampled by sampler 68, an equalizer 70 
compensates for the distortion introduced by LPF 64 and the local loop, and then a 
decision device or decoder 72 estimates the transmitted level by selecting the level 
that is closest to the received point. From the level the decoder 72 determines the 

25 equivalence class and then recovers the six information bits by performing an 
inverse mapping function. 

The operation of receiver 66 is essentially unchanged as compared to the 
receiver described in the co-pending applications referred to above. The only 
difference is that the receiver now needs to consider a larger set of possible levels 

30 and the inverse mapping involves the determination of the equivalence class. 

Equalizer 70 compensates for the linear distortion introduced by the LPF 64 and the 
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local loop 65, as described in the co-pending applications. For example, when a 
linear equalizer is used, the output of the equalizer can be represented as follows: 



rk = yk + nk (1 ) 

where nk is the total noise plus distortion present at the output of the equalizer. 
Decoder 72 then selects the levels yk nearest to rk as the decision, determines its 
equivalence class, and then recovers the six information bits by an inverse map. 

If the equalizer includes a maximum-likelihood sequence estimator (e.g., the 
Viterbi equalizer), then the received signal can be represented in the form 

rk = X yk-j fj + nk, (2) 

and this time, the decoder selects the closest sequence {yk} using a Viterbi decoder. 
For each estimated symbol yk, the decoder determines its equivalence class and 
then finds the six information bits via an inverse map. 

Encoder 56, FIG. 6, includes MAP 74 which is a look-up table containing for 
each possible combination of the six-bit groups of data received from bit collector 54, 
FIG. 5, levels representing each equivalence class i, where i is an integer between 0 
and 63. Each level, two in this example, y(i,1) and y(i,2) is provided to level selector 
76 where a decision is made as to which level, yk, is to be transmitted. 

This decision is made as follows. First, encoder 56 keeps track of the running 
digital sum (RDS) of the transmitted levels, yk, by feeding back the output of level 
selector 76 to function block 78. From the previously transmitted levels, yk, function 
block 78 calculates the weighted RDS, zk = -(1-b)RDS, where 0 s= b < 1 is weighting 
factor. Because of D/A nonlinearities, the exact values of the yk levels may not be 
known in encoder 56; however, this should not have a significant effect. It is possible 
to determine the error and send this information back to encoder 56 to make these 
calculations more accurate. 

Given the group of six bits to be transmitted, level selector 76 selects as the 
level yk from the equivalence class {y(i,1), y(i,2)} the level closest to the weighted 
RDS. It can be seen that when the RDS is positive, zk will be negative and vice 
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versa. This enables the encoder to choose a level, yk, from each equivalence class 
such that when its value is added to the RDS it will bring it closer to zero than the 
other levels in the equivalence class. After selecting the level yk the octet which 
represents the level yk is determined by octet converter 80 and transmitted over the 
5 digital network. The value of the transmitted octet can be obtained from a look-up 
table. 

The variable b is a weighting factor that controls the trade-off between the 
sharpness of the spectral null and the average energy of the transmitted signal. Our 
analysis has shown that when the number of levels is sufficiently larger than the 

10 number of equivalence classes, the sequence yk will have a spectrum which can be 
approximated by the filter response h(D) = (1-D)/(1-bD). Clearly, when b = 0, we find 
that h(D) = 1-D, which is the well-known Class I Partial Response with a sinusoidal 
spectral shape having a null at DC. On the other hand, as b approaches 1 , the 
spectrum becomes flat across much of the band except for a very sharp spectral null 

15 at DC. It can be seen that for b = 0, the average energy of yk will be twice as large 
as in the case of a flat spectral shape. As b approaches 1 , however, the average 
energy increase will disappear. In some applications, it may be desirable to keep the 
constellation expansion, measured by the ratio of the number of levels to the number 
of equivalence classes. 

20 It will be apparent to those skilled in the art that the invention can be used with 

constellations of any number of levels, and with any smaller number of equivalence 
classes. 

The present invention may be more broadly utilized to spectrally shape, as 
desired, the analog signals output from the |n-law to linear converter at the central 

25 office. The example described above is a specific case of using this invention to 
reduce the energy of the transmitted signal around DC, but the principals of this 
invention used in that example can be generalized to spectrally shape signals in 
numerous ways, for example, to pre-equalize the signals. 

A generic version of the encoder of this invention, encoder 56a, is shown in 

30 FIG. 7. The only difference between this general case and the special case of a 
spectral null described above is how the sequence or spectral function zk is 
generated. Let h{D) be a monic, causal impulse response of a filter representing the 

9 
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desired spectral shape, where D is a delay operator. Suppose we represent the 
sequences {yk} and {zk} using D-transform notation as y(D) and z(D), respectively. 
Then, the sequence z(D) can be represented as 



z(D) = (1 - 1/h(D)) y(D). (3) 

A close examination of this equation reveals that at a given time k, zk only 
depends on past values of yk, and therefore can be determined recursively. Thus, 
for each six bit group, encoder 56a determines which level from the associated 
equivalence class is nearest in value to zk and selects that level. The octet 
representing that level is then transmitted. Again, our analysis shows that for 
sufficiently large number of levels the sequence {yk} transmitted by the central office 
60 will have a spectrum closely approximating the spectrum of the filter with 
response h(D). 

The technique described here can also be used in conjunction with a more 
complex scheme for mapping the information bits to equivalence classes. For 
example, it can be used in conjunction with shell mapping, a mapping technique 
used in the V.34 high-speed modem specification. 

The examples described above are for an uncoded system. However, the 
principals can be easily applied to a coded system, for example a trellis coded 
system. The only difference in this case is that the equivalence classes are further 
partitioned into subsets, which are used to construct the trellis code. 

For example, when a one-dimensional trellis code based on a 4-way set 
partition is utilized together with the same 64-level signal constellation to send 5 bits 
per symbol, the equivalence classes are partitioned into subsets as follows: ai , bi , 
C1 , di , a2, b2, C2, d2, ...a n , b n , c n ,d n . In the example described above, the 64 
equivalence classes would be partitioned into four subsets each containing sixteen 
equivalence classes. The output of a rate-1/2 convolutional encoder, e.g. two of the 
six bits in a group, then determines the subset, and the remaining four "uncoded" bits 
select the specific equivalence class within the subset. The actual level from the 
chosen equivalence class in the chosen subset is selected as described above. The 
operation of the encoder is otherwise unchanged. 

10 
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Of course, when trellis coding is utilized, the receiver will use a decoder to 
select the most likely sequence. The trellis decoder may also be an equalizer, jointly 
decoding the trellis code and equalizing for intersymbol interference. 

It may also be possible to use the present invention to enable detection of 
loss of frame synchronization in a receiver. This can be accomplished by 
infrequently, but periodically violating the rule for selecting the signal point in a given 
equivalence class, where the period is chosen to be an integer multiple of the 
desired framing. A loss of frame synchronization, can be detected in the receiver by 
monitoring such rule violations. The receiver can also reacquire frame 
synchronization or may simply request a synchronization pattern (training sequence) 
from the transmitter. 

PCM Upstream Precodina 

There is shown in FIG. 8, a typical PCM communication system 100. System 
100 includes analog PCM modem 102 connected to a telephone company central 
office (CO) 104 over a local analog loop or channel 103. There is also included a 
digital network 106 which is interconnected to CO 104 and to digital PCM modem 
1 08. With this system, PCM data may be transmitted both in the downstream 
direction (i.e., from digital PCM modem 108 to analog PCM modem 102) and in the 
upstream direction (i.e., from analog PCM modem 102 to digital PCM modem 108). 
This type of bi-directional PCM communication system is described in US Application 
Serial No. 08/724,491, entitled Hybrid Digital/Analog Communication Device, which 
is assigned to the assignee of the present invention and which is incorporated herein 
in its entirety by reference. 

In the above section a technique for PCM downstream spectral shaping or 
precoding of data signals is described. In this section there is described a 
precoding technique for PCM upstream precoding of data signals. 

In FIG. 9 there is shown in block diagram 1 10, an example of PCM upstream 
transmission in accordance with this invention. In block diagram 1 10 there is 
included analog PCM modem 112 interconnected to analog channel 113. Analog 
PCM modem 112 includes transmitter 120 having a precoder 122, prefilter 124 and a 
digital to analog converter (D/A) 126. Precoder 122 receives digital data u(n) and 

n 
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outputs precoded digital data signal x(n). The precoded digital data signal is filtered 
by prefilter 124 to form signal z(n) which is provided to D/A 126. D/A 126 converts 
the filtered signal z(n) to analog form and transmits analog signal, z(t), over analog 
channel 1 13, having a channel characteristic, c(t). 

The analog channel modifies the transmitted signal z(t) to form signal y(t). 
The signal y(t) then encounters downstream PCM echo, echo(t) 128, that is added to 
y(t), producing signal r(t). Signal r(t) is received by u.-law (A-law in some countries 
outside of the US) quantizer 130 in central office (CO) 1 14 and is quantized 
according to the in-law. See International Telecommunications Union, 
Recommendation G.71 1, Pulse Code Modulation (PCM) of Voice Frequencies, 1972. 

The quantized octets (digital values), q(n), are transmitted over digital network 
1 1 6 at a frequency of 8kHz where they may be affected by various digital 
impairments, as discussed below. The possibly affected octets, v(n), are received by 
digital PCM modem 1 18 which ideally decodes the octets, v(n), into their 
corresponding constellation points, y(t), from which the original digital data, u(n), can 
be recovered. The decoding of v(n) is described in co-pending application entitled 
System, Device and Method for PCM Upstream Transmission Utilizing an Optimized 
Transmit Constellation, CX097028, which is assigned to the assignee of the present 
invention and which is incorporated herein in its entirety by reference. 

Before data can be transmitted upstream, the clock (f,) of D/A 126 in analog 
PCM modem 112 must be synchronized to the clock (f 2 ) of CO 1 14. This can be 
achieved by learning the clock from the downstream PCM signal (not shown) and 
synchronizing the clocks using the technique proposed in US Patent No. 5,199,046, 
entitled First and Second Digital Rate Converter Synchronization Device and 
Method, incorporated herein by reference in its entirety. Once the clocks are 
synchronized, PCM upstream block diagram 110, FIG. 9, can be represented as 
equivalent discrete time block diagram 1 10', FIG. 10, with like components being 
represented by the same reference numbers containing a prime ('). In block diagram 
1 1 0' it is assumed that f, = f 2 ; however, it must be noted that f, does not have to be 
equal to f 2 as long as the two clocks are synchronized. When f, is equal to f 2 ,n is the 
time index for 8kHz samples, since the clock (f 2 ) of CO 24 is fixed at that frequency. 
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An example where f t does not equal f 2 is depicted in FIG 1 1 . Equivalent 
discrete time block diagram 1 10a\ FIG. 11, is the same as equivalent discrete time 
block diagram 1 10\ FIG. 10, except that there is a 2X up-sampler 123a' in 
transmitter 120a' and a 2X down-sampler 129a' to account for the fact that i,=2f 2 . 
The variables "m" and "n" are the time indexes for 16kHz and 8kHz samples, 
respectively. 

Precoder 122' and prefilter 124 1 , according to this invention, are designed to 
transmit signal z(n) over analog channel 1 1 3 such that predetermined constellation 
points, y(n), corresponding to digital data u(n) are produced at the input of \i-\aw 
quantizer 130' (in combination with an echo component, echo(n), if present). In other 
words, the input of |n-law quantizer 130* is y(n) + e(n) in the presence of echo(n) and 
just y(n) in the absence of echo(n). 

Using the PCM upstream preceding technique described below, or another 
precoding technique, it is difficult for digital PCM modem 1 18' to accurately decode 
u(n) from v(n) in the presence of echo, quantization and digital impairments without a 
properly designed transmit constellation of points, y(n). It is described in co-pending 
application CX097028 how to design the transmit constellation for y(n) to enable y(n) 
(and eventually u(n)) from v(n)) to be decoded in the presence of echo, quantization 
and digital impairments with minimized error probability. 

As described in co-pending application CX097028, for a given connection, 
depending on the line conditions, a transmit constellation for each robbed bit 
signaling (RBS) time slot is selected. As an example, transmit constellation 140 is 
depicted in FIG. 12. This constellation includes ten constellation points, y 0 -y 9 , 
ranging in value from -39 to 39. It should be noted that the constellation points, y(n), 
are not necessarily G.71 1 jx-law levels. 

The constellation points y(n) correspond to digital data to be transmitted, u(n). 
In other words, each constellation point represents a group of data bits and the 
number of data bits represented by each constellation point depends on the number 
of points in the constellation (and the number of equivalence classes which are 
described below). The more points in the constellation, the more bits of data that 
can be represented. As shown in FIG. 12, digital data u(n) is divided into four groups 
of bits 0,1 ,2 and 3, corresponding to 00, 01 , 10 and 1 1, for example. Thus, in this 
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example each constellation point transmitted represents two bits and since the 
constellation points are transmitted at 8k/sec, the data rate is 16kbps. It must be 
understood that this is a simplified example and data may be mapped into u(n) using 
any mapping schemes that can map bits into equivalence classes, such as shell 
mapping or modulus conversion. 

According to this invention, the constellation points are grouped into 
equivalence classes. An equivalence class is a set of typically two or more 
constellation points which represent the same group of bits or digital data to be 
transmitted, u(n). With constellation 140, it is shown that constellation points y 0 (-60), 
y<(-6), and y 8 (45) form the equivalence class for u(n)=0. Constellation points y,(-45), 
y 5 (6), and y 9 (60) form the equivalence class for u(n)=1 and constellation points y 2 (- 
31), and y 6 (18) form the equivalence class for u(n)=2. Finally, constellation points 
y 3 (-18), and y 7 (31) form the equivalence class for u(n)=3. 

Equivalence class selection is generally accomplished as follows. The 
constellation, with M points, is indexed as y 0 , y,....^ in ascending (or descending) 
order. Assuming u(n) has U values, e.g. U=4 as in the above example, then the 
equivalence class for u(n)=u contains all the y k 's where k modulo U is u. For 
example, in FIG. 1 1 , the equivalence class for u(n)=0 is y 0 , y u , y^, where U=4. Note 
that each equivalence class is not required to have the same number of constellation 
points. 

The number of supporting data levels for u(n) should be chosen to satisfy the 
following two conditions: 1) The expansion ratio, which is defined as the ratio 
between the number of constellation points for y(n) and the number of supporting 
data levels for u(n), i.e., M/U; and 2) TX power constraints. 

The expansion ratio should be large enough to guarantee stable operation. 
The size of the expansion ratio will depends on the channel characteristics. In voice 
band modem applications, there is at least one spectral null at f=0. Therefore, we 
should have an expansion ratio of M/U > 2 to make the system stable. In practice, to 
guarantee the stability, the quality of the channel is determined from the channel 
response, c(n), and the minimum expansion ratio is set accordingly. For example, 
we can use C(f=4kHz), the frequency response of the channel at 4kHz (with respect 
to other frequencies like 2kHz), as the quality of the channel and depending on that 
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quality we set the minimum expansion ratio. If the C(J = 4kHz) * C{J = 2kHz) , then 



we set M/U> 2.0 . As the C(/ = 4kHz) gets smaller and smaller, the expansion 
ratio must be increased. 

As described below, precoder 122' selects the appropriate constellation point, 
y k , from the equivalence class for the data, u(n), to be transmitted and determines a 
value for x(n) that will produce the selected constellation point at the input to u.-law 
quantizer 1 30' . 

The precoding scheme, i.e., the design of precoder 122' and prefilter 124*, are 
now described as follows. From the characteristics of analog channel 113', c(n), 
n=0,1,...N c -1, determined by digital PCM modem 118', as described in co-pending 
application entitled Device and Method for Detecting PCM Upstream Digital 
Impairments in a Communication Network, CX097029, which is assigned to the 
assignee of the present invention and which is incorporated herein in its entirety by 
reference, an optimal target response p(n), n=0,1,...N p -1, and corresponding prefilter 
g(n), n=-A,-A+1,..., -A+N B -1 (where A is the decision delay), as shown in FIG. 10, are 
determined. This problem is similar to determining the optimal feedforward and 
feedback filters for a decision feedback equalizer (DFE). The prefilter corresponds to 
feedforward filter of DFE and the target response corresponds to feedback filter of 
DFE. See, N. Al-Dhahir, et al, "Efficient Computation of the Delay Optimized Finite 
Length MMSE-DFE", IEEE Transactions On Signal Processing, vol. 44, no. 5, May 
1996, pp.1 288-1 292. Preferably, the target response p(n) and the filter g(n) will be 
determined in the analog modem, but they can be determined in the digital modem 
and transmitted to the analog modem. 

The prefilter g(n), n=-A,-A+1 -A+N B -1 , and the target response p(n), 

n=0,1,...,N p -1, (where p(0)=1) can be derived given c(n) by minimizing the cost 
function £ as follows: 



The first term ensures small intersymbol Interference (ISI), i.e., the receiver of 
digital PCM modem 118' receives what precoder 122' tried to encode, and the 




(4) 
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second term enforces the transmit (TX) power to stay finite and small. The term a is 
a constant term which should be chosen depending on the application. The larger a 
is the lower the TX power will be, but at the expense of ISI. A smaller a will give less 
ISI at the expense of TX power. Therefore a should be chosen depending on what is 
desired for ISI and TX power for a given application. As an example, a can be 
chosen to be the signal to noise ratio (SNR) of the system, which is <j\ /E(x 2 ) or 

SNR normalized by channel energy, i.e., SNR^c f . For E(x 2 ), we can use -9dBm 
which is the power constraint for upstream transmission. This minimization problem 
is the same as DFE tab initialization problem. The term o„ 2 can be determined as 
described in co-pending application CX097028. 

The initially determined p(n) and g(n) can always be used if the analog 
channel c(n) is time invariant. However, in practice, c(n) is time variant, though it is 
very slowly changing. Therefore, some kind of adaptation scheme is necessary. One 
way to do it is to monitor performance and retrain if the performance goes bad, i.e., 
re-estimating c(n) in the digital modem 1 18' and sending a new c(n) back to analog 
modem 112* to recalculate g(n) and p(n). Another way is to feedback the analog 
channel error signal, error(n), as described in co-pending application CX097029, 
from digital modem 1 18* to analog modem 112' through downstream data 
transmission and use that error signal to adapt p(n) and g(n). 

Once the target response p(n) is determined precoder 122' can be 
implemented. As explained above, we can send data u(n) by transmitting x(n) such 
as to produce at the input to quantizer 130", FIG. 10, a constellation point y(n) which 
is one of the points in the equivalence class of u(n). Which constellation point from 
the equivalence class of u(n) to use to represent u(n) is usually selected to minimize 
the TX power of transmitter 120'. The TX power of transmitter 120' is the power of 
z(n) (or some other metric). In practice, since it is hard to minimize the power of z(n), 
the power of x(n) is minimized instead, which is a close approximation of minimizing 
z(n). 

The following is a known relationship among x(n), y(n) and p(n): 



y(n) = p(n) * x(n) 
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where "*" represents convolution. That relationship can be expressed as follows: 

y(n) = p(0)x(n) + p(1 )x(n-1 ) + ... p(N p )x(n-N p ) (6) 
Since p(0) is designed to equal to 1 , then equation (6) can be simplified as follows: 

Mn)=y(n)-^j^i)x(ri-i) . (7) 
«'=i 

And, since p(n) and the past values of x(n) are known, the appropriate y(n), 
among the constellation points of the equivalence class of a given u(n), may be 
selected to minimize x 2 (n) in order to minimize the TX power of transmitter 120*. 

Or, lookahead (i.e., decision delay) can be introduced to choose y(n). That is, 
y(n-A) can be chosen from the set of equivalence classes for u(n-A) to minimize 
I x(n-A) | 2 + | x(n-A+1 ) 1 2 + ...+ I x(n) I a , where: 

x(n -j) = y(n -j) - ^p(i)x(n -j-i) (8) 
i=i 

where j=0,1 .... A and where y(n-j) is chosen from the set of equivalence classes of 
u(n-j) 0=0,1,... A-1). - 

Precoder 122' may be implemented according to this invention as depicted in 
FIG. 13. Precoder 122' includes a mapping device 150 which receives the incoming 
digital data u(n) from a digital data source and, depending on the number of bits that 
can be transmitted with each constellation point, determines for each group of bits 
the equivalence class associated with the group of bits. Mapping device 1 50 outputs 
the constellation points, y k , forming the equivalence class to TX signal/constellation 
point selector 152 which selects the constellation point, y k .from the equivalence 
class and determines the transmit signal x(n) based on the input from calculation 
device 154. 

Filter device 154 receives the transmit signal x(n) and calculates the 
summation term (or running filter sum (RFS)) of equation (7) above. Based on the 
value of the RFS, TX signal/constellation point selector 152 selects the constellation 
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point in the equivalence class that will cause x(n) in equation (7) to be closest in 
value to zero and calculates the value of x(n) from the calculated RFS and the 
selected constellation point. The calculated transmit signal x(n) is then provided to 
prefilter 124' where x(n) is filtered to form signal z(n) which is transmitted over analog 
channel 113', FIG. 10. 

In order to limit the TX power of transmitter 120', FIG. 10, to keep it within the 
FCC regulations, the equivalence classes for u(n) must be designed accordingly. 
With a constellation having a predetermined number of constellation points, If we 
want to send more data, then more groups of data, u(n), and hence equivalence 
classes for u(n) will be required. As a result, the constellation points will be further 
away and will require more transmit power. This is because y(n) is chosen as 
described below according to equation (7) to minimize x 2 (n). Therefore, if the 
constellation points in the equivalence classes are spaced further apart, it is more 
likely that x*(n) will be larger. Thus, to reduce the TX power, we can make the 
equivalence class of u(n) closer at the expense of rate. This is depicted in FIGS. 14A 
and 14B. 

In FIGS. 14A and 14B, both constellations 156, FIG. 14A, and 158, FIG. 14B, 
have the same number of constellation points; however, constellation 1 56 has only 
three equivalence classes u(n)= 0,1 and 2 while constellation 158 has five 
equivalence classes u(n)= 0,1,2,3 and 4. Using constellation 158 will require more 
TX power than constellation 156, but it will be capable of transmitting at a higher data 
rate. 

The approximate TX power (the power of z(n) ) can be calculated as follows 
when U is the number of points desired to support u(n): 

j U-l 

P z =\g(nf — X dist\u{n)= /)/12 (9) 
U i=o 

where lg(n)l 2 is the energy of prefilter and dist(u(n)=i) is the minimum distance 
between the points in the equivalence points. For example, in FIG. 12 dist(u(n)=0) = 
I -6-(-60) | = 54. Several values of U should be tried to find out the one which 
satisfies the power constraints. Note also that this should be done for each time slot. 
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The transmit constellation selection and equivalence class selection according 
to this invention may be summarized as follows: 

1 ) Obtain digital impairments, calculate noise variance, ct„ 2 , and echo variance, 
a e 2 , as described in co-pending application CX097028; 

2) From c e 2 , a n 2 , and the digital impairments, choose the proper constellation for 
y(n) for each time slot, also as described in co-pending application CX097028 ; and 

3) For each time slot, find the number of points that can be supported for u(n) 
while satisfying the TX power constraints and the minimum expansion ratio to 
guarantee stable operation. From this U the constellation for y(n), and the 
equivalence classes for u(n) can be determined. 

The above precoding technique which utilizes a one dimensional constellation 
can be expanded to multi-dimensional constellations by expanding the definition of 
the equivalence class of u(n). The following references describe various 
downstream precoding techniques using multi-dimensional constellations: Eyuboglu, 
Vedat; "Generalized Spectral Shaping for PCM Modems," Telecommunications 
Industry Association, TR30.1 Meeting, Norcross, Georgia, 9-11 April 1997, pages 1- 
5; Eyuboglu, Vedat; "Convolutional Spectral Shaping," Telecommunications Industry 
Association, TR30.1 Meeting, Norcross, Georgia, 9-11 April 1997; Eyuboglu, Vedat; 
"More on Convolutional Spectral Shaping," ITU Telecommunications Standardization 
Sector 009, V.pcm Rapporteur Meeting, La Jolla, CA, 5-7 May 1997; Eyuboglu, 
Vedat; "Draft Text for Convolutional Spectral Shaping," ITU-T SG 1 6 Q23 
Rapporteur's Meeting, September 2-1 1 , 1997, Sun River, Oregon; Eyuboglu, Vedat; 
"A Comparison of CSS and Maximum Inversion," Telecommunications Industry 
Association, TR30.1 Meeting on PCM Modems, Galveston, Texas, 14-16 October 
1997; and Eyuboglu, Vedat; "Draft Text for Convolutional Spectral Shaping," 
Telecommunications Industry Association, TR30.1 Meeting Galveston, Texas, 14-16 
October 1997. 

Moreover, the example described above is for an uncoded system. However, 
the principals can be easily applied to a coded system, for example a trellis coded 
system. The only difference in this case is that the equivalence classes are further 
partitioned into subsets, which are used to construct the trellis code. 



19 



WO 99/34565 PCT/US98/24368 

Generalized PCM Precodina 

The above described PCM upstream precoding technique (i.e. from analog 
PCM modem 1 12', FIG. 10, to digital PCM modem 118, may be applied to an analog 
PCM modem to analog PCM modem connection as depicted in FIG. 15. System 160 
includes analog PCM modem 162 connected to CO 166 over analog loop or channel 
164. CO 166 is interconnected to digital network 168. Similarly analog PCM modem 
174 is connected to CO 170 over analog loop or channel 172. And, CO 170 is 
connected to digital network 168. 

Block diagram 180, FIG. 16, depicts an analog PCM modem to analog PCM 
modem connection according to this invention. In block diagram 180 there is 
included analog PCM modem 182 interconnected to analog channel 184. Analog 
PCM modem 182 includes transmitter 200 having a precoder 202, prefilter 204 and a 
digital to analog converter (D/A) 206. Precoder 202 receives digital data u(n) and 
outputs precoded digital data x(n). The precoded digital data is filtered by prefilter 
204 to form signal z(n) which is provided to D/A 206. D/A 206 converts the filtered 
signal z(n) to analog form and transmits analog signal, z(t), over analog channel 
184, having a channel characteristic, c(t). 

The analog channel modifies the transmitted signal z(t) to form signal y(t). 
The signal y(t) then encounters PCM echo, echo(t) 208, that is added to y(t), 
producing signal r(t). Signal r(t) is received by fi-law (A-law in some countries 
outside of the US) quantizer 210 in central office (CO) 186 and is quantized 
according to the u-law. See International Telecommunications Union, 
Recommendation G.711, Pulse Code Modulation (PCM) of Voice Frequencies, 1972. 

The quantized octets (digital values), q(n), are transmitted over digital network 
1 88 at a frequency of 8kHz where they may be affected by various digital 
impairments, as discussed below. The possibly affected octets, v(n), are received by 
CO 1 90 and the octets, v(n), are converted by n-law D/A 212 into analog levels for 
transmission over analog channel 1 92. The levels are received by analog PCM 
modem 194 which converts the levels to data u(n). 

Once the clocks f1 to f2 of D/A 206 and D/A 210 are synchronized, block 
diagram 180 can be modeled as discrete time block diagram 180', FIG. 17. Analog 
PCM modem should do the equalization to get v(n) from g(n) in the same way as a 
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downstream PCM modem works as is known in the art. Then, from v(n), a PCM 
upstream decoding algorithm to decode y(n), i.e. u(n), is undertaken. 

The above only describes transmission from analog PCM modem 1 82* to 
analog PCM 1 94'; however, transmission in the other direction is accomplished in the 
same manner. The above described PCM upstream precoding technique (i.e. from 
analog PCM modem 1 12', FIG. 10, to digital PCM modem 118,) can be applied 
directly to an analog PCM modem to analog PCM modem connection as depicted in 
FIGS. 15-17. 

It should be noted that this invention may be embodied in software and/or 
firmware which may be stored on a computer useable medium, such as a computer 
disk or memory chip. The invention may also take the form of a computer data 
signal embodied in a carrier wave, such as when the invention is embodied in 
software/firmware which is electrically transmitted, for example, over the Internet. 

The present invention may be embodied in other specific forms without 
departing from the spirit or essential characteristics. The described embodiments 
are to be considered in all respects only as illustrative and not restrictive. The scope 
of the invention is, therefore, indicated by the appended claims rather than by the 
foregoing description. All changes which come within the meaning and range within 
the equivalency of the claims are to be embraced within their scope. 
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1 . A transmitter for precoding a sequence of analog levels transmitted over an 
analog channel to a quantization device, wherein the analog channel modifies the 
transmitted analog levels, the transmitter comprising: 

a mapping device for mapping data bits to be transmitted to a sequence of 
equivalence classes, wherein each equivalence class contains one or more 
constellation points; and 

a constellation point selector interconnected to the mapping device which 
selects a constellation point in each equivalence class to represent the data bits to 
be transmitted and which transmits a level that produces the selected constellation 
point at an input to the quantization device. 

2. The transmitter of claim 1 further including a filter device, operably coupled to 
the constellation point selector, which receives at its input previously transmitted 
levels and provides its output to the constellation point selector. 

3. The transmitter of claim 2 wherein the constellation point selector selects the 
constellation point from each equivalence class based on the output of the filter 
device. 

4. The transmitter of claim 3 further including a prefilter, having a predefined filter 
response, g(n), for filtering the level transmitted by the constellation point selector. 

5. The transmitter of claim 4 wherein the response of the filter device is: 

p(i)x(n - i) 

where p(i) is a target response and x(n-i) represents the previously transmitted 
levels. 

6. The transmitter of claim 5 wherein the target response, p(n), and the prefilter 
response, g(n), are derived from the predetermined response, c(n), of the analog 
channel. 
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7. The transmitter of claim 5 wherein the constellation point selector transmits 
the levels, x(n), according to the following function: 



5 where y(n) are the constellation points. 

8. The transmitter of claim 7 wherein the constellation point selector selects the 
constellation point in each equivalence class which minimizes the transmit power of 
the transmitter by selecting the constellation point, y(n), which produces the smallest 
10 value for x(n). 
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9. A method for transmitting a precoded sequence of analog levels over an 
analog channel to a quantization device, wherein the analog channel modifies the 
transmitted analog levels, the method comprising: 

mapping data bits to be transmitted to a sequence of equivalence classes, 
wherein each equivalence class contains one or more constellation points; and 

selecting a constellation point in each equivalence class to represent the data 
bits to be transmitted; and 

transmitting a level that produces the selected constellation point at an input 
to the quantization device. 

10. The transmitter of claim 9 wherein the step of selecting a constellation point 
includes filtering the previously selected constellation points with a filter device and 
selecting the constellation points based on the output of the filter device. 

1 1 . The method of claim 1 0 further including filtering the level transmitted with a 
prefilter having a predefined filter response, g(n). 

1 2. The method of claim 1 1 wherein the response of the filter device is: 

Z,p(i)x(n-i) 

where p(i) is a target response and x(n-i) represents the previously transmitted 
levels. 

13. The method of claim 12 wherein the target response, p(n), and the prefilter 
response, g(n), are derived from the predetermined response, c(n), of the analog 
channel. 

1 4. The method of claim 1 2 wherein step of transmitting includes transmitting the 
levels, x(n), according to the following function: 

^n)=y(n)-2jii)x(n~i) 
i=i 

where y(n) are the constellation points. 
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15. The method of claim 14 wherein the step of selecting includes selecting the 
constellation point in each equivalence class which minimizes the transmit power of 
the transmitter by selecting the constellation point, y(n), which produces the smallest 
value for x(n). 
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16. A computer useable medium having computer readable program code means 
embodied therein for transmitting a precoded sequence of analog levels over an 
analog channel to a quantization device, wherein the analog channel modifies the 
transmitted analog levels, the method comprising: 

computer readable program code means for mapping data bits to be 
transmitted to a sequence of equivalence classes, wherein each equivalence class 
contains one or more constellation points; and 

computer readable program code means for selecting a constellation point in 
each equivalence class to represent the data bits to be transmitted; and 

computer readable program code means for transmitting a level that produces 
the selected constellation point at an input to the quantization device. 

17. The computer useable medium of claim 16 wherein the computer readable 
program code means for selecting a constellation point includes computer readable 
program code means for filtering the previously selected constellation points with a 
filter device and selecting the constellation points based on the output of the filter 
device. 



1 8. The computer useable medium of claim 1 7 further including computer 
readable program code means for filtering the level transmitted with a prefilter 
having a predefined filter response, g(n). 

1 9. The computer useable medium of claim 18 wherein the response of the filter 
device is: 

p(i)x(n - f ) 

where p(i) is a target response and x(n-i) represents the previously transmitted 
levels. 

20. The computer useable medium of claim 19 further including computer 
readable program code means for deriving the target response, p(n), and the prefilter 
response, g(n), from the predetermined response, c(n), of the analog channel. 
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21 . The computer useable medium of claim 19 wherein the computer readable 
program code means for transmitting includes computer readable program code 
means for transmitting the levels, x(n), according to the following function: 

n p 

xtri) = y(n)-^^i)x(n-i) 
i=i 

5 where y(n) are the constellation points. 

22. The computer useable medium of claim 21 wherein the computer readable 
program code means for selecting includes computer readable program code means 
for selecting the constellation point in each equivalence class which minimizes the 

10 transmit power of the transmitter by selecting the constellation point, y(n), which 
produces the smallest value for x(n). 
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23. A computer data signal embodied in a carrier wave, wherein embodied in the 
computer data signal are computer readable program code means for transmitting a 
precoded sequence of analog levels over an analog channel to a quantization 
device, wherein the analog channel modifies the transmitted analog levels, the 
method comprising: 

computer readable program code means for mapping data bits to be 
transmitted to a sequence of equivalence classes, wherein each equivalence class 
contains one or more constellation points; and 

computer readable program code means for selecting a constellation point in 
each equivalence class to represent the data bits to be transmitted; and 

computer readable program code means for transmitting a level that produces 
the selected constellation point at an input to the quantization device. 

24. The computer data signal of claim 23 wherein the computer readable program 
code means for selecting a constellation point includes computer readable program 
code means for filtering the previously selected constellation points with a filter 
device and selecting the constellation points based on the output of the filter device. 

25. The computer data signal of claim 24 further including computer readable 
program code means for filtering the level transmitted with a prefilter having a 
predefined filter response, g(n). 

26. The computer data signal of claim 25 wherein the response of the filter device 
is: 



where p(i) is a target response and x(n-i) represents the previously transmitted 
levels. 

27. The computer data signal of claim 26 further including computer readable 
program code means for deriving the target response, p(n), and the prefilter 
response, g(n), from the predetermined response, c(n), of the analog channel. 
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28. The computer data signal of claim 26 wherein the computer readable program 
code means for transmitting includes computer readable program code means for 
transmitting the levels, x(n), according to the following function: 

x(ri) = y(n)-2jKt)rtn-i) 
where y(n) are the constellation points. 

29. The computer data signal of claim 28 wherein the computer readable program 
code means for selecting includes computer readable program code means for 
selecting the constellation point in each equivalence class which minimizes the 
transmit power of the transmitter by selecting the constellation point, y(n), which 
produces the smallest value for x(n). 



29 



WO 99/34565 PCTAJS98/24368 

30. In an analog pulse code modulation (PCM) modem adapted for upstream 
PCM data transmission to a digital PCM modem, a transmitter for precoding a 
sequence of analog levels transmitted over an analog channel to a quantization 
device, wherein the analog channel modifies the transmitted analog levels, the 
transmitter comprising: 

a mapping device for mapping data bits to be transmitted to a sequence of 
equivalence classes, wherein each equivalence class contains one or more 
constellation points; and 

a constellation point selector interconnected to the mapping device which 
selects a constellation point in each equivalence class to represent the data bits to 
be transmitted and which transmits an analog level that produces the selected 
constellation point at an input to the quantization device. 
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31 . In an analog pulse code modulation (PCM) modem adapted for PCM data 
transmission to another analog PCM modem, a transmitter for precoding a sequence 
of analog levels transmitted over an analog channel to a quantization device, 
wherein the analog channel modifies the transmitted analog levels, the transmitter 
comprising: 

a mapping device for mapping data bits to be transmitted to a sequence of 
equivalence classes, wherein each equivalence class contains one or more 
constellation points; and 

a constellation point selector interconnected to the mapping device which 
selects a constellation point in each equivalence class to represent the data bits to 
be transmitted and which transmits an analog level that produces the selected 
constellation point at an input to the quantization device. 
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